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By Seyed G. Saddoughi 
1. Motivation and background 

This is a report on the continuation of our experimental investigations (Sad- 
doughi 1993; Saddoughi & Veeravalli 1993) of the hypothesis of local isotropy in 
shear flows. This hypothesis, which states that at sufficiently high Reynolds num- 
bers the small-scale structures of turbulent motions are independent of large-scale 
structures and mean deformations (Kolmogorov 1941), has been used in theoretical 
studies of turbulence and computational methods such as large-eddy simulation. 
Since Kolmogorov proposed his theory, there have been many experiments, con- 
ducted in wakes, jets, mixing layers, a tidal channel, and atmospheric and labora- 
tory boundary layers, in which attempts have been made to verify - or refute - 
the local-isotropy hypothesis. However, a review of the literature over the last five 
decades indicated that, despite all these experiments in shear flows, there was no 
consensus in the scientific community regarding this hypothesis, and, therefore, it 
seemed worthwhile to undertake a fresh experimental investigation into this ques- 
tion. 


1.1 Plane boundary layer 

In our previous reports, we presented hot-wire measurements of the velocity 
fluctuations in the test-section-ceiling boundary layer of the 80- by 120-foot Full- 
Scale Aerodynamics Facility at NASA Ames Research Center, the world’s largest 
wind tunnel. At our measurement location, the boundary- layer thickness, <5, was 
about 1.1m, and the maximum Reynolds numbers based on momentum thickness, 
Re, and on Taylor microscale, R\, were approximately 370,000 and 1,450 respec- 
tively. These were the largest ever attained in laboratory boundary- layer flows. 
The boundary layer developed over a rough surface, but the Reynolds-stress pro- 
files agreed with canonical data sufficiently well for our purposes. Spectral and 
structure-function relations for isotropic turbulence were used to test the local- 
isotropy hypothesis, and our results established the condition under which local 
isotropy can be expected. 

Here we use a Cartesian co-ordinate system Xj = (x, y, z) with jr-axis along 
the flow direction, y-axis normal to the solid surface, and ^-axis in the span- 
wise direction. The respective mean-velocity components in these directions are 
Ui = (U,V,W), and the fluctuating components are u t = (u,u,u>). Overbars de- 
note time averages. 

Our plane boundary- layer data showed that, to within the accuracy of measure- 
ment, the shear-stress co-spectral density £^ 12 (^ 1 ), which is the most sensitive indi- 
cator of local isotropy, fell to zero at a wavenumber about a decade larger than that 
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at which the energy spectra first followed -5/3 power laws. At the highest Reynolds 
number, £ 12 (^ 1 ) vanished about one decade before the start of the dissipation range, 
and it remained zero in the dissipation range. 

We found that the lower-wavenumber limit of locally-isotropic behavior of the 
shear-stress co-spectra is given by k\yjej S 3 ~ 10 where S is the mean shear, 
dU/dy, and e is the average turbulent energy dissipation rate per unit mass. Our 
investigation also indicated that for energy spectra this limit could be relaxed to 
kiy/ffs* « 3; this is Corrsin’s (1958) criterion, with the numerical value ob- 
tained from our data. The existence of an isotropic inertial range requires that 
this wavenumber be much less than the wavenumber at the onset of viscous effects, 
jbjrj < 1 (?? is the Kolmogorov length scale), so that the combined condition (Corrsin 
1958 and Uberoi 1957) is Sy/vje C 1. 

Among other detailed results, it was observed that in the dissipation range, the 
energy spectra had a simple exponential decay (Kraichnan 1959) with an exponent 
prefactor close to the value 0 = 5.2 obtained in direct numerical simulations at low 
Reynolds number. Plots of compensated spectra, k s / 3 ,E(ki), proved to be a very 
sensitive test in the inertial subrange. The Kolmogorov constants obtained from 
the one-dimensional data at high Reynolds numbers satisfied the isotropic relations 
for the spectra and the second-order structure functions, and the constant for the 
three-dimensional spectrum, C, was estimated to be 1.5 ± 0.1 (Monin &: Yaglom 
1975). Spectral “bumps” between the -5/3 inertial range and the dissipative range 
were observed on all the compensated energy spectra. The shear-stress co-spectra 
rolled-off with a -7/3 power law and scaled linearly with 5 (Lumley 1967). 

In summary, our results confirmed the local-isotropy hypothesis for simple shear 
layers, and it was shown that one decade of iner tial subrange with truly negligible 
shear-stress co-spectral density requires S^/vJe not more than about 0.01 (for a 
shear layer with turbulent kinetic energy production « dissipation, this implies a 
microscale Reynolds number of about 1500). 

1.2 Distorted boundary layer 

The effects of extra mean strain rates on the large-scale structure of shear flows 
(Bradshaw 1973) have been investigated extensively. The unanswered question for 
us was “will our criteria for the existence of local isotropy hold for complex non- 
equilibrium flows”. Therefore, experiments to address this question were designed 
for the 80- by 120-foot wind tunnel. One possible experiment was to study the 
plane-of- symmetry flow in front of an obstacle placed vertically in a fully- developed 
two-dimensional turbulent boundary layer, e.g. a circular cylinder placed with its 
axis perpendicular to the plate. There have been a number of experimental inves- 
tigations dealing with the large-scale structural changes that occur in this kind of 
flow (Johnston 1960; Hornung & Joubert 1963; Belik 1973; Mehta 1984, Agui & 
Andreopoulos 1990; Devenport & Simpson 1990; to name a few). In this type of 
boundary layer, the pressure rises strongly as the obstacle is approached and in the 
plane of symmetry of the flow the boundary layer is also influenced by the effects 
of lateral divergence (Saddoughi Joubert 1991). Hence, in addition to the basic 
shear, dU/dy, the extra strain rates involved in the flow are dU/dx, dV/dy and 
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dW/fc To obtain the desired effects, the size of the cylinder should be at least of 

sn' K T SS ° f thC b ° Undary layer ' To conduct such an experiment 

m the 80 by 120 wind tunnel, a very large cylinder had to be fixed to the ceiling 

ot the tunnel. This presented considerable construction difficulties. Here some of 

the results taken in front of this cylinder at the highest speed of the tunnel will be 
discussed. 


2. Accomplishments 

2.1 Apparatus and measurement techniques 
The only possible way for attaching an obstacle to the ceiling of the 80- by 120- 
foot wind tunnel was to use one of the existing light ports for this purpose. Since the 
diameter of our test cylinder had to be larger than the clear opening of a typical 
light port, we had to use two concentric cylinders. As shown in Plate 1(a), the 
main cylinder is a ready-made light- weight polyethylene tank (wall thickness = I", 
diameter D = 4' and height L = 6') and the inner (second) cylinder is a A" thick 
8 diameter aluminum (6061-T651) tube, which extends for about 6' into the attic 
through the light port. Steel rods are bolted to the aluminum cylinder. A thick 

high-density polyethylene plate and a thick aluminum plate are respectively 
bolted to the inside and outside of the bottom of the polyethylene tank The 
effective wall thickness at the bottom of the tank is §". Another thick high- 
density polyethylene plate is bolted to the top of the tank, and the space between 
the aluminum cylinder and the tank is filled with high-density Polyurethane foam. 
This provided us with a fairly light-weight (500 lb) and solid cylinder. 

From the attic above the test-section ceiling, this whole unit was pulled up 
through the light port and was attached to the attic structure. Plate 1(b) shows 
the test cylinder in place as viewed from inside the test section. The NASA safety 
requirements were satisfied, and the cylinder and its attachments were designed for 
maximum possible aerodynamic load and dynamic loading due to seismic activity 
with an overall factor of safety of 5 on yield strength. 

Another light port upstream of the cylinder was used for traversing the probe 
through the boundary layer. Therefore, our measurement location was fixed at 
x/D « 0.85 with respect to the front of the cylinder. The measurement strategy 

^^ men ! a r, tl0n ’ and P rocedure were ^ similar to those explained by Saddoughi 
(1993) and Saddoughi & Veeravalli (1993), and details will not be repeated here. 

2.2 Results and discussion 

The measurements to be discussed here correspond to the maximum reference 
ve ocity of the tunnel, U ref = 51.25 m/s. The normalized profiles of the longitudi- 
nal mean velocity, U/Urtj , for the present distorted boundary layer are compared 
with the profiles obtained for both the high-speed and low-speed cases of the plane 
boundary layer in figure 1, where y is the distance from the wall. Also shown in 
this figure are the least-squares polynomial fit to each of the data sets (solid lines) 
which have been used to obtain the mean-flow integral parameters for each experi- 
ment. The shape of the velocity profile for the distorted boundary layer is typical 
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Plate 1. Test cylinder, (a) Outer (main) cylinder (polyethylene tank: wall thick- 
ness = i", diameter D = 4', and height L = 6') and the inner (support) cylinder 
(aluminum 6061-T651 tube: thick, 8" diameter and 12' height) (b) Test ^ cylin- 

der attached to the ceiling of the 80' X 120' wind tunnel, as viewed from inside the 

test section. 
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y(m) 

FIGURE 1. Normalized longitudinal mean-velocity profiles measured in plane 
and distorted boundary layers, o , U re f * 50 m/s and □ , U ref « 10 m/s plane 
boundary layer; ® , U re f ss 51 m/s distorted boundary layer. The solid lines are 
the least-square polynomial fit to each data set. 

of the adverse-pressure-gradient flows: reduction in dU/dy. The boundary-layer 
thickness, <5 (the point where U/U e = 0.995) has increased to approximately 1250 
mm in the distorted boundary layer. Here the shape factor H as 1.85, and in the 
freestream the pressure coefficient C p « 0.23. 

Figure 2 shows the normalized profiles of the normal velocity component, V/U re f, 
and the spanwise component, W/U ref , which were measured by the X-wires in 
UV- and UW-mode respectively. It can be seen that, as expected, in the plane of 
symmetry of the flow the crossflow, W, is approximately equal to zero. A least- 
squares polynomial fit to the V profile was used to obtain the values of dV/dy. 

The magnitudes of the extra strain rate due to the streamline divergence, QiW j dz, 
influencing the plane of symmetry of the flow can be obtained from (9W/ dz) — 
U(d/3/dz) (see e.g. Saddoughi & Joubert 1991), where (3 is the flow yaw angle mea- 
sured at different spanwise locations z. The profiles of j3 measured by a yaw-meter 
probe for three spanwise locations (z/D = -0.21,0,0.21) through the boundary 
layer are shown in figure 3. The profiles are typical of three-dimensional bound- 
ary layers: larger flow yaw angles near the wall than the freestream. Finally, the 
continuity equation was used to obtain the <917 jdx values. 

The profiles of the Reynolds normal stresses («f/^ e/ ,u|/f/^ / ,u|/t/2 ), and 
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FIGURE 2. Normalized vertical and spanwise mean- velocity profiles measured in 
the plane of symmetry of the distorted boundary layer. 0 , V /U rt j^ W/U T ef • 



FIGURE 3. Flow yaw-angle profiles measured in the distorted boundary layer at 
different spanwise locations. 
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y (m) 

FIGURE 4. Profiles of Reynolds stresses measured in plane and distorted bound- 
ary layers: (a) ^f/[/ re / 2 , (b) u| l/U ref 2 , (c) u 2 3 /U r ef 2 , (d) uTu^/Uref 2 - For key to 
symbols see figure 1. 
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Figure 5. The structure parameter, aj = —uiu 2 /q 2 . The solid line is the 
canonical smooth-wall value (Townsend 1976). For key to symbols see figure 1. 


the shear stress, —uTu 2 /U 2 e j, for the distorted and plane boundary layers are com- 
pared in figure 4. The profiles for the distorted boundary layer appear to be quite 
different from the plane flow case. The peaks of u| and the shear stress, — tri'u 2 , 
profiles have moved away from the wall to y ss 300 mm, and in the outer part 
of the layer the values of all the Reynolds stresses have increased. The changes 
in the large-scale structure of turbulence can be seen in figure 5, which shows the 
profiles of Townsend’s structure parameter, aj = (— uiiz 2 /g 2 ), where q 2 (= TTpuJ) is 
twice the turbulent kinetic energy per unit mass. The large drop in the values of 
this parameter in the inner part of the boundary layer, and their recovery to the 
canonical values in the outer part of the layer are apparently due to the effects of 
adverse pressure gradients (see Bradshaw 1967). 

The spectral measurements of the three components of the velocity made at 
y — 100 mm, 300 mm, and 500 mm are analyzed here. The location y — 300mm 
was chosen because, as shown earlier, the peak of the Reynolds stresses occurred 
at this position in the layer. For these three locations, the values of the extra- 
mean-strain rates, ( dU/dx)/(dU/dy ), ( dV / dy) / (dU / dy ), and ( dW/dz)/(dU /dy ) 
are all larger than 0.1, which according to Bradshaw (1973) should produce large 
non-linear effects on the large-scale structures of the boundary layer. 

In figure 6, we examine the II( = -^b tJ b Jt ) and III( = f bijb jk b ki ) invariants map 
of the Reynolds stress anisotropy tensor (Lumley &: Newman 1977; Lee & Reynolds 
1985), bij = ( UiUj/q 2 ) — Sijj 3, where 6ij = 1 or 0 for i = j or i ^ j respectively. 
The data points shown in this figure correspond to those positions in both the 
boundary layers where spectral measurements were taken. It appears that for both 
the plane and distorted layers, the values are close to the axisymmetric expansion 
limit. However, note that there is a fairly large increase for the distorted boundary 
layer, particularly at y = 100 mm, which is represented by the highest point in each 
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FIGURE 6. Anisotropy invariant map of Reynolds stress. Open and solid symbols 
axe for the plane and distorted boundary layers respectively. The maximum value 
in each case corresponds to measurements at y — 100 mm. 


data set. 

As an example of the spectral measurements, figure 7 shows Kolmogorov’s uni- 
versal scaling of the one-dimensional longitudinal power spectra at the inner-layer 
position of the distorted boundary layer and the mid-layer position of the plane 
boundary layer, compared with a compilation of previous experimental work taken 
from Chapman (1979) with later additions. The Reynolds numbers for the present 
distorted boundary layer have increased substantially, and a maximum R\ ~ 1960 
has been obtained for the inner-layer position. 

The intensity (or rapidity) of a mean strain rate can be measured in terms of 
s = y/sijS t] /2, (Lee & Reynolds 1985). To be consistent with our earlier definition, 
we shall use S = 2s as the equivalent mean strain rate. To test the local-isotropy 
hypothesis in the distorted layer, log-linear plots of the correlation- coefficient spec- 
tra, R 12 ( h) = -E l2 (h)l are plotted versus the non-dimensional 
wavenumber Aq s/e/S 3 in figure 8. In isotropic flow the shear-stress co-spectrum, 
Ei 2 {k i), which satisfies / 0 °° £ 12 (^ 1 ) dki = — u\u 2 , is equal to zero. This indicates 
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Figure 7. Kolmogorov’s universal scaling for one-dimensional longitudinal power 
spectra. The present mid-layer and inner-layer spectra for the respective plane and 
distorted boundary layers are compared with data from other experiments. This 
compilation is from Chapman (1979), with later additions. 
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FIGURE 8. Correlation- coefficient spectra obtained at different locati ons in the 
distorted boundary layer for U re f = 51 m/s: wavenumber scaled with y/e/S 3 . (a) 
y = 500 mm, R\ « 1960; (b) y = 300 mm, R\ « 1950; (c) y = 100 mm, i?A « 1750. 
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FIGURE 9. Ratios of calculated to measured transverse spectra at different loca- 
tions and freestream velocities in t he pla ne ( ) and distorted ( ) boundary 

layers: wavenumber scaled with yjej S 3 . 

that for local isotropy the correlation- coefficient spectrum should fall to zero at 
high wavenumbers. For all the thr ee m easurement positions in the boundary, the 
Corrsin-Uberoi condition S * H S yfv I £ 1 was satisfied. It can seen from figure 8 

that at al l the measurement locations J?i 2 (fci) reaches the isotropic value of zero at 

\J Tj 5 3 w 10, which is the same limit found for the plane boundary layer for the 
onset of local isotropy. 

If the motion is isotropic, the transverse spectra £*22(^1 ) and E§${k \ ) are uniquely 
determined by the longitudinal spectrum E n (k x ) (e.g. Batchelor 1953): £ 22 (fc 1) = 
£33(^1) = |(1 - ~) E n(k\)' The transverse spectra, E$ c (k\) and 

can be calculated from the measured longitudinal spectrum, F^ eas (iti), using the 
above equation. An anisotropy measure may be defined as E™J c (ki )/E™* a9 (k \ ), 
where <y = 2 or 3 corresponds to u 2 or U3 respectively. These anisotropy measures 
should be equal to 1.0 in an isotropic flow. We have used 9^-order, least-squares 
polynomial log-log fits to the data to calculate these measures, which are shown 
in figure 9 scaled using the length scale y/e/S 3 . The data for the plane boundary 
layer are also shown in this figure. It is obvious that the uncertainty in estimat- 
ing the 5 values for the distorted boundary layer is larger than the plane case. 
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FIGURE 10(a). For caption see next page. 

However, it appears that a reasonably good collapse is obtained and local isotropy 
of energy spectra is achieved (±10%) in the inertial subrange for non-dimensional 

wavenumbers k\ yj ej S 3 >3. . 

For both plane and distorted boundar y laye rs, the collapse of the longitudinal 
and tr ans verse spectra achieved using y/eJS* and y/T[S as length and velocity 
scales, respectively, is shown in figure 10. As to be expected, these scales will not 
collapse the low and high-wavenumber ranges of the spectra. Also marked on this 
figure is the wavenumber corresponding to the start of Rn{ki) ~ 0, which is about 
one decade higher than the start of the -5/3 law on the energy spectra. This 
plot clearly demonstrates the fact that a —5/3 law does not necessarily imply local 
isotropy. We also note that in the high-wavenumber range the extent of the -5/3 
law does not increase with R x , but it is a function of (e/S 2 u), which is the Reynolds 
number based on the above length and velocity scales. However, the accuracy of 
this observation will be examined further after completing the upcoming low-speed 
measurements . 
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FIGURE 10 . Spectra measured at different locations and freestream velocities in the 
plane and distorted boundary layers non-dimensionalized using y/e/S* and yfzJS 
as length and velocity scales respectively, (a) (b) £22(^1); (c) £33(^1). 
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